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I.  SUMMARY  OF  WORK  AND  DESCRIPTION  OF  WORK  ACCOMPLISHED 


This  project  has  involved  the  closely  coordinated  theoretical  and  experimental 
investigation  of  the  collisional  relaxation  of  the  molecular  free  radicals  OH  and  NO  by  open- 
shell  atoms.  These  species  are  generated  in  high  vibrational  and  rotational  levels  in  the 
mesosphere  and  lower  thermosphere  through  highly  exothermic  reactions,  e.g.  N(2D)  +  02  and  H 
+  03.  The  atmospheric  signatures  of  these  and  other  molecules  are  strongly  influenced  by  the 
subsequent  redistribution  of  this  energy.  While  there  has  been  considerable  study  of  the 
relaxation  of  OH  and  NO  by  stable  molecules,  our  knowledge  of  their  relaxation  by  open-shell 
atoms,  notably  0(3/>),  prior  to  this  work  was  quite  limited. 

Efforts  under  this  project  have  mainly  focused  on  the  study  of  collisional  relaxation  of 
OH  and  NO  by  open-shell  atoms,  in  particular  0(3P).  The  sections  below  summarize  the 
accomplishments  of  this  project,  as  well  as  describe  work  initiated  with  earlier  support  under  the 
HEDM  program.  For  administrative  reasons,  the  collaboration  between  Alexander  and 
Dagdigian  was  supported  by  two  separate  grants,  one  to  the  University  of  Maryland  and  one  to 
The  Johns  Hopkins  University. 

A.  Theoretical  Investigation  of  Collisions  of  a  Free  Radical  with  an  Open-Shell  Species 

A  statistical  model,  originally  developed  by  Rackham  and  Manolopoulos1'2  to  describe 
chemical  reactions  proceeding  through  a  strongly  bound  collision  complex,  has  been  applied  to 
the  collisional  relaxation  of  a  free  radical  by  an  open-shell  species.  In  such  a  collision,  at  least 
one  of  the  potential  energy  surfaces  possesses  a  deep  well  corresponding  to  the  stable  molecule 
which  can  be  formed  by  recombination  of  the  open-shell  species.  The  inelastic  scattering 
dynamics  is  treated  by  coupled-states  capture  theory  with  inclusion  of  all  the  relevant  potential 
energy  surfaces  and  the  couplings  between  them.  Specifically,  we  have  shown  how  to  introduce 
the  coupling  between  the  nuclear  and  electronic  degrees  of  freedom  of  a  molecule  in  a  2n 
electronic  state.. 

We  have  used  this  approach  to  investigate  the  inelastic  scattering  of  OH  by  H.3  Since  the 
O('D)  +  H2  reaction  is  exothermic,  at  energies  below  the  O^D)  +  H2  channel  threshold  capture 
(from  OH  +  H)  into  the  H20  complex  can  lead  only  to  release  of  OH  in  a  different  vibration- 
rotation  state.  The  overall  inelastic  scattering  can  be  divided  into  a  direct  component  due  to 
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coupling  between,  and  the  anisotropy  of,  the  repulsive  potentials  which  correlate  with  OH  +  H 
and  a  capture  component,  due  to  the  capture-release  mechanism  discussed  above. 

Our  calculations3  show  that  rotationally  inelastic  scattering  [OH(/'j  +  H  — »OH(/')  +  H]  is 
dominated  by  direct  scattering  and  is  similar  in  magnitude  and  dependence  on  final  quantum 
number  of  the  rotationally  inelastic  scattering  of  OH  by  He  or  H2.  However,  for  vibrational 
relaxation  the  capture  mechanism  dominates  and  results  in  vibrationally  inelastic  cross  sections 
which  are  not  insubstantial.  This  confirms  an  earlier  conjecture  by  Smith4  of  complex-mediated 
vibrational  relaxation.  In  addition,  we  predict  the  same  final-state  OH  A'/ A"  propensities  seen  in 
our  theoretical  simulations  of  the  O  +  H2  reaction,5  as  well  as  in  the  earlier  experimental  work  of 
Wiesenfeld,  Butler,  Gericke,  and  their  co-workers.6  7 

Potential  energy  surfaces  (PESs)  for  the  interaction  of  OH  with  0(3P)  atoms  have  also 
been  computed.8  Both  the  lowest  2 A'  and  2 A"  potential  energy  surfaces  have  a  barrierless  access 
to  a  well.  In  contrast  to  the  O('D)  +  H2  reaction,  here  the  H  +  02  channel  is  lower  in  energy  than 
the  OH  +  O  asymptote.  Consequently,  the  complex-mediated  process  OH(v=l)  +  O  — ►  OHO+  — * 
products  pathway  leads  to  removal  of  the  vibrationally  excited  level  primarily  by  production  of 
H  +  02  products.  Rate  constants  have  been  computed  using  statistical,  coupled-states 
calculations2  for  collisions  on  the  2 A'  and  2 A"  potential  energy  surfaces,  separately.  Overall  good 
agreement  was  seen  between  the  calculated  thermal  vibrational  relaxation  rate  constant  at  300K 
and  the  value  measured  in  Dagdigians’s  laboratory,9  which  is  discussed  in  the  next  section. 
However,  for  the  OH(v=0)  +  O  — *•  H  +  02  reaction  a  comparison  between  the  statistical, 
coupled-states  integral  cross  sections  and  the  more  exact,  full  reactive  scattering  calculations  of 
Honvault  and  Guo."1  indicates  that  the  statistical  method  overestimates  the  reactive  cross  section 
by  at  least  a  factor  of  50%.  We  are  currently  continuing  our  theoretical  investigation  of  this 
reaction. 

B.  Experimental  Investigation  of  Collisional  Relaxation  of  Vibrationally  Excited  OH  and 
NO  by  Open-Shell  Atoms 

Knowledge  of  rate  constants  for  collisional  relaxation  of  the  internal  degrees  of  freedom 
of  the  OH  and  NO  radicals  are  important  for  an  understanding  of  the  IR  atmospheric  emissions 
of  these  species.  While  there  has  been  considerable  study  of  relaxation  rates  by  stable,  closed- 
shell  collision  partners,  only  a  few  investigations  of  collisions  with  open-shall  species  have  been 
carried  out.  We  have  carried  out  experimental  studies  of  the  collisional  vibrational  relaxation  of 
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OH(v=l)  by  0(3P)9  and  N(4,S)"  atoms  and  have  investigated  the  collisional  vibrational  relaxation 
of  NO  by  oxygen  atoms. 

Hydroxyl  radicals  in  the  ground  v=0  vibrational  level  can  react  with  both  0(3P)  and 
N(4S)  atoms,  to  yield  02  and  NO  molecular  products,  respectively,  and  OH(v=l)  can  be  removed 
by  both  reaction  and  vibrational  relaxation.  For  both  atomic  species,  collisional  removal  of 
OH(v=l)  is  mediated  by  formation  and  decay  of  a  transient  complex  (H02  and  HNO, 
respectively).  Our  experiment  employs  a  discharge-flow  apparatus,  in  which  O (3P)  or  N(4S) 
atoms  are  prepared  by  microwave  discharge  in  an  02/Ar  or  N2/Ar  mixture,  respectively.  A 
photolytic  precursor  is  added  upstream  of  a  laser-induced  fluorescence  detection  port,  and  OH  is 
generated  by  193  nm  photolysis  of  a  suitable  precursor  (>95%  H202).  We  had  considered  the 
use  of  the  commonly  employed  OH  precursor  HNO3,  but  193  nm  photolysis  produces  some 
O^D)12  and  complicates  the  relaxation  kinetics  by  its  reaction  with  the  precursor  to  yield  OH  in 
several  excited  vibrational  levels.  OH(v=0  and  1)  concentrations  were  monitored  by  laser- 
induced  fluorescence  as  a  function  of  time  after  the  photolysis  laser  pulse.  Photolysis  of  H202 
yields  mainly  OH(v=0)  and  -1%  OH(v=l),13  and  possible  cascade  contributions  from  higher  OH 
vibrational  levels  are  avoided. 

In  these  experiments,  the  concentration  of  the  atoms  was  determined  from  the  decay  rate 
of  the  OH(v=0)  concentration  when  the  microwave  discharge  is  on.  This  is  a  reasonable 
approach  for  the  measurement  of  the  atom  concentrations  since  the  rate  constants  for  the 
OH(v=0)  +  0(3P)  — >  H  +  02  and  OH(v=0)  +  N(4S)  — >  H  +  NO  reactions  are  known  from 
previous  work.14’ 15  In  effect,  we  are  comparing  the  collisional  removal  of  OH(v=0  and  1)  under 
the  same  conditions  in  order  to  determine  the  enhancement  of  the  removal  rate  with  vibrational 
excitation.  In  extracting  the  total  rate  constants  for  removal  of  OH(v=l)  in  collisions  with  the 
atoms,  we  must  also  take  into  account  vibrational  relaxation  of  OH(v=l)  by  H202  and  other 
species.  In  the  case  of  collisions  with  0(3P),  we  found  that  the  concentration  of  the  H202 
precursor  in  the  flow,  as  gauged  by  the  OH(v=0  and  1)  laser-induced  fluorescence  signals,  was 
reduced  by  a  factor  of  ~3  when  the  microwave  discharge  is  turned  on.  This  results  from  a  slow 
chain  reaction  involving  0(3P),  H02,  OH,  and  H  radicals  which  consumes  H202.  We  simulated 
this  decomposition  by  kinetic  modeling  in  order  to  estimate  species  concentrations  in  the  low 
with  the  microwave  discharge  on  and  hence  0(3P)  atoms  present.  In  the  case  of  OH-O 
collisions,  we  determined  that  the  ratio  of  the  room-temperature  total  removal  rate  constant  for 
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OH(v=l)  to  that  for  OH(v=0)  was  1.1  ±0.1.  Thus,  the  removal  rate  for  OH(v=l)  is  found  to  be 
slightly  larger  than  for  OH(v=0),  but  the  difference  is  within  the  experimental  uncertainty. 

Schatz  and  co-workers16  carried  out  an  extensive  ab  initio  study  of  the  PESs  connecting 
the  OH  +  N  and  H  +  NO  channels,  and  the  higher-energy  O  +  NH  asymptote.  The  OH  +  N 
reagents  can  react  on  a  3A"  PES  to  form  a  HON  complex,  which  can  isomerizes  to  form  HNO 
and  then  dissociate  to  H  +  NO  products.  Unfortunately,  they  did  not  explore  the  dependence  of 
the  3  A "  PES  upon  the  O-H  separation.  Following  on  our  study  of  0H(v=l)-0  collisions,  we 
studied  the  effect  of  vibrational  excitation  on  the  total  removal  rate  of  OH(v=l)  by  N(4S) 
atoms.11  In  this  system,  the  OH  concentration  was  unaffected  by  the  presence  or  absence  of 
N(4S)  atoms.  For  this  system,  we  found  that  the  ratio  of  the  room-temperature  total  removal  rate 
constant  for  OH(v=l)  to  that  for  OH(v=0)  was  1.6±0.4.  Hence,  there  is  a  significant 
enhancement  in  the  removal  rate  with  vibrational  excitation  in  this  system. 

We  also  investigated  the  vibrational  relaxation  of  OH(v=l)  by  H  atoms,  which  were 
prepared  by  a  microwave  discharge  of  H2.  In  this  case,  the  OH  concentrations  were  reduced  by  a 
factor  of  ~10  when  the  discharge  was  turned  because  of  a  chain  reaction  of  the  precursor  H202 
with  H  atoms,  much  greater  than  was  the  case  for  0(3P)  atoms.  The  greater  efficiency  of 
removal  of  H202  by  H  atoms  is  consistent  with  our  kinetic  modeling  of  the  destruction  of  H202. 
Because  of  the  significant  destruction  of  H202  by  H  atoms,  a  reliable  rate  constant  for  the 
collisional  removal  of  OH(v=l)  by  H  atoms  could  not  be  determined. 

We  also  carried  out  an  experimental  study  of  collisional  vibrational  relaxation  of  NO  by 
0(3P)  atoms.  The  NO(v)  concentration  is  followed  as  a  function  of  time  after  the  193  nm 
photolysis  of  r-butyl  nitrite  precursor  by  laser-induced  fluorescence  detection  in  the  NO  A  -  X 
band  system.  In  this  experiment,  the  oxygen  atom  concentration  is  determined  by  cavity  ring- 
down  spectroscopy  on  the  forbidden  JD2  <—  3P2  line,  as  described  in  detail  in  Sec.  D.  With  this 
concentration  measurement,  the  pseudo  first-order  NO(v)  decay  rate  can  be  converted  to  a 
bimolecular  rate  constant,  after  correction  for  other  (small)  loss  processes.  A  preliminary  value 
of  (2.86±0.61)  x  10-11  cm3  molecule-1  s-1  was  obtained  for  the  relaxation  of  NO(v=2)  by  O(  V) 
atoms.  This  value  falls  in  the  lower  range  of  recent  measurements  of  the  room-temperature 
vibrational  relaxation  rate  constants  of  NO(v=l  and  2)  by  oxygen  atoms. 
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C.  Investigation  of  Electronic  Quenching  of  OH(A2E+)  in  Collisions  with  Molecular 
Hydrogen 

In  a  collaborative  project  with  the  Lester  group  at  the  University  of  Pennsylvania,  we 
carried  out  theoretical  calculations  to  help  understand  the  experimental  investigation  of  the 
inelastic  scattering  channel  resulting  from  electronic  quenching  of  OH  A2Z+  by  molecular 
hydrogen  (also  supported  by  AFOSR  under  grant  FA9550-04- 1-0068).  The  experiments  utilized 
a  pump-probe  scheme  to  determine  the  OH  X  2H  population  distribution  following  collisional 
quenching  in  a  pulsed  supersonic  expansion.  17  18  The  experiments  revealed  a  high  degree  of 
rotational  excitation  (A")  of  the  quenched  OH  X  2Tl  products  observed  in  v"=l  and  2,  as  well  as 
a  pronounced  propensity  for  quenching  into  the  FI(A')  A-doublet  level.  To  supplement  these 
experiments,  we  carried  out  extensive  multi-reference,  configuration-interaction  calculations 
aimed  at  exploring  the  topology  of  the  relevant  potential  energy  surfaces.  Electronic  quenching 
of  OH  A  2Z+  by  H2  proceeds  through  conical  intersections  between  two  potentials  of  A'  reflection 
symmetry  (in  planar  geometry)  that  correlate  with  the  electronically  excited  A2Z+  and  ground  X 
2U  states  of  OH.  The  conical  intersections  occur  in  high-symmetry  geometries,  in  which  the  O- 
side  of  OH  points  toward  H2.  Corroborating  earlier  work  of  Hoffman  and  Yarkony,19  our 
calculations  revealed  a  steep  gradient  away  from  the  OH-H2  conical  intersection  as  a  function  of 
both  the  OH  orientation  and  interfragment  distance.  The  former  is  responsible  for  the  high 
degree  of  rotational  excitation  observed  for  the  quenched  OH  X  2Yl  products.  In  addition, 
because  a  third  PES  of  A"  reflection  symmetry,  correlating  with  OH  X2!!  +  H2,  does  not  have 
this  pronounced  angular  gradient,  coplanar  passage  through  the  conical  intersection  should  favor 
the  observed  A-doublet  propensity.  In  addition,  we  predict  that,  as  the  quenched  products 
separate,  adiabatic  relaxation  of  the  H2  stretch  will  be  incomplete,  leading  to  substantial 
vibrational  excitation  of  the  H2  products  subsequent  to  the  electronic  quenching  of  OH. 

D.  Determination  of  0(3P)  Atom  Concentrations  by  Cavity  Ring-Down  Spectroscopy  on 
the  Forbidden  *D2  <—  3/>2  Line 

Oxygen  atoms  are  important  in  many  environments,  including  combustion  media, 
plasmas,  and  the  atmosphere.  The  availability  of  a  non-intrusive  spectroscopic  method  for  the 
determination  of  0(3P)  atom  concentrations  is  highly  desirable.  We  desire  such  a  method  in  our 
collisional  relaxation  experiments.  The  measurement  of  rate  constants  for  collisional  energy 
transfer  in  free  radicals  induced  by  oxygen  atoms  entail  measurement  of  the  decay  rates  of  the 
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free  radical  concentration  under  pseudo  first-order  conditions,  i.e.  with  the  oxygen  atoms  in 
excess.  The  desired  bimolecular  rate  constants  can  then  be  calculated  by  dividing  the  decay  rates 
by  the  oxygen  atom  concentration. 

Several  spectroscopic  methods  have  previously  been  employed  to  determine  O f3/5)  atom 
concentrations.  The  strong  3 s  3S  <—  3P  resonance  transition  near  130  nm  has  been  extensively 
employed  in  discharge-flow  experiments20  but  is  limited  to  low  concentrations  and  cannot  be 
applied  to  flames  etc.  Two-photon  fluorescence  excitation21  near  226  nm  to  the  3 p  3P  state  and 
detection  of  near  IR  3 p  3P  — »  35  3S  emission  circumvents  problems  with  VUV  spectroscopy; 
however,  this  approach  requires  consideration  of  electronic  quenching,  energy  transfer, 
stimulated  emission,  and  ionization,  as  well  as  requiring  independent  calibration.22  24 

We  have  shown  that  the  use  of  cavity-ring  down  spectroscopy  (CRDS)25  26  on  the 
forbidden  *D2  3P2  l*ne  at  630  nm  can  be  employed  for  the  non-intrusive  spectroscopic 

determination  of  Of3/5)  atom  concentrations.27  This  magnetic  dipole  transition  has  a 
recommended28  spontaneous  emission  probability  of  5.63  x  10-3  s_1,  which  implies  a  peak 
absorption  cross  section  of  5.71  x  10-22  cm2  for  a  room-temperature  Doppler-broadened  line 
profile,  or  an  integrated  absorption  cross  section  of  2.93  x  10-23  cm  molecule-1.  In  CRDS, 
absorption  is  detected  by  a  increase  in  the  rate  of  decay  of  light  injected  into  an  optical  cavity 
containing  the  sample. 

We  have  applied  CRDS  to  detect  0(3/>)  atoms  in  a  discharge-flow  apparatus.  Oxygen 
atoms  were  generated  by  passing  a  N20/Ar  mixture  or  pure  02  through  a  2450  MHz  microwave 
discharge.  The  effluent  from  the  discharge  passed  through  a  104  cm  long  optical  cavity,  83  cm 
of  which  contained  Of3/5)  atoms.  The  CRD  mirrors  have  99.99%  reflectivity,  which  led  to  a 
empty-cavity  photon  decay  lifetime  of  ~35  /ns.  In  order  to  take  account  of  the  fact  that  our  laser 
linewidth  (0.15  cm-1  FWHM)  was  greater  than  the  O  atom  Doppler  width  (0.05  cm-1  FWHM), 
integrated  O  atom  absorbances  were  determined  by  scanning  over  the  line.  We  found  that  Of3/5) 
atom  concentrations  in  the  10 14  cm-3  range  could  be  measured  to  an  accuracy  of  ~4  x  1013 
cm-3. 

There  is  an  uncertainty  of  only  8%  in  the  transition  probability  for  the  *D2  3P2  line,28 

and  hence  a  similar  uncertainty  in  the  absorption  cross  section.  We  checked  the  accuracy  of  our 
CRDS  measurements  by  comparison  concentrations  determined  by  this  technique  and  by  the 
established  method  of  titration  with  N02  and  found  good  agreement  between  the  two  methods. 
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We  are  currently  employing  this  method  to  determine  0(3P)  atom  concentrations  in  our  ongoing 
study  of  the  vibrational  relaxation  of  NO(v)  by  0(3P)  atoms,  as  mentioned  above. 

E.  Simulations  of  A1  Atoms  Doped  in  Liquid  Helium 

We  described  for  publication  the  results  of  a  path- integral,  Monte-Carlo  simulation  of  the 
spectrum  of  a  single  A1  atom  doped  into  liquid  helium  clusters.29  For  this  system  high-resolution 
spectra  at  0.38K  have  been  recorded  by  Lehmann,  Scoles,  and  coworkers,  with  AFOSR  HEDM 
support.30  The  helium  environment  is  unique  because  the  interaction  between  helium  atoms  is 
very  weak.  Because  of  the  large  zero-point  motion,  helium  clusters  are  known  to  be  the  only 
clusters  that  remain  liquid  under  all  conditions  of  formation.31  Although  the  spectra  of  alkali 
atoms  in  helium  droplets  have  been  investigated  theoretically  by  density  functional  methods  (Cs) 
and  path-integral  Monte-Carlo  simulations  (Li,  Na,  K),31'32  our  work  is  the  first  attempt  to 
perform  similar  simulations  for  group  Ilia  atoms. 

Following  earlier  work  in  our  group,33  36  we  used  ab  initio  Al-He  pair  potentials  for  the 
ground  and  lowest  excited  Al-He  electronic  states,  within  a  Balling  and  Wright  pairwise 
Hamiltonian  model,37  to  describe  the  interaction  between  the  open-shell  A1  atom  and  an  arbitrary 
number  of  He  atoms.  By  describing  this  interaction  as  a  sum  of  pairwise  6x6  Hamiltonian 
matrices,  rather  than  the  sum  of  pairwise  scalar  potentials,  we  can  treat  accurately  the  electronic 
anisotropy  of  the  A1  atom.  A  similar  approach  is  used  to  describe  the  interaction  of  A1  in  its 
excited  3 d  and  4 p  states  with  multiple  He  atoms. 

In  our  simulation  of  the  arrangement  of  He  atoms  around  a  single  A1  atom  doped  in  a  He 
cluster,  we  used  a  multilevel  Monte-Carlo  scheme  to  minimize  the  computational  effort.  This 
involves  first  performing  large  coarse  displacements  of  the  system  in  configuration  space,  and 
then  refining  only  if  the  initial  displacement  is  accepted,  an  adaptation  of  the  bisection  algorithm 
of  Ceperly  and  co-workers.3s  The  simulation  of  the  absorption  spectra  was  done  using  a 
semiclassical  Franck-Condon  approximation.39,40 

Our  calculations  reveal  that  the  doping  of  the  A1  3 p  electron  strongly  influences  the  He 
packing.  With  inclusion  of  tail  corrections  for  the  ground  and  excited  states  potentials,  our 
calculated  3z/< — 3/?  spectrum  agrees  reasonably  well  with  the  experimental  spectrum  reported  by 
Lehmann,  Scoles,  and  co-workers.30  The  blue  shift  of  the  calculated  spectrum  associated  with 
the  4s  <—  3 p  transition  of  solvated  A1  is  about  25  nm  (~2000  cm-1)  larger  than  seen  in 
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experiments  on  A1  embedded  in  bulk  liquid  He.41  We  also  predict  that  the  spectrum  associated 
with  the  4 p  <—  3 p  transition  will  be  blue  shifted  by  ~7000  cm-1  (nearly  1  eV). 

F.  Recombination  of  A1  Atoms  in  Solid  Hydrogen 

The  embedding  of  light  impurities  in  solid  hydrogen  has  received  considerable 
experimental  and  theoretical  attention  in  recent  years,  partly  because  of  the  potential 
technological  applications  as  high  energy  density  materials.42  In  collaboration  with  the  Voth 
group  at  Utah,  we  have  published  several  earlier  investigations  of  a  single  B35  or  Al36  atom 
embedded  in  solid  para- H2.  We  found  that  in  both  cases  the  dependence  of  the  atom-/?H2 
potential  on  the  orientation  of  the  singly-filled  p  electron  induces  some  distortion  of  the  nearest- 
neighbor  ligand  shell  to  allow  the  /7-orbital  to  orient  itself  to  minimize  the  repulsive  I  and 
maximize  the  attractive  n  interactions. 

For  open-shell  atoms  such  as  B  and  Al,  one  important  question  remains  unanswered: 

when  doped  with  multiple  atomic  impurities,  is  the  system  stable  with  respect  to  recombination, 

a  highly  exothermic  process?  To  study  more  than  one  doped  open-shell  atom,  we  must  consider 

the  interaction  between  impurity  atoms,  which  depends  on  the  relative  orientation  of  the  p- 

electrons  of  both  dopants.  In  an  earlier  article  we  developed  a  theoretical  framework  for  the 

2 

accurate  description  of  the  interaction  of  two  P  atoms  in  the  presence  of  multiple  spherical 
ligands.43 

With  this  potential,  we  developed44  a  simple  approach  to  investigate  the  likelihood  of 
recombination:  We  place  the  two  Al  atoms  at  several  initial  substitution  sites,  similar  to  those 
chosen  in  our  earlier  paper.43  We  then  allow  the  separation  vector  for  the  two  Al  atoms  to 
decrease.  At  each  reduced  value  of  this  separation  coordinate,  which  we  keep  fixed,  we  then 
perform  path-integral,  Monte-Carlo  calculations  to  sample  the  motion  of  the  p\\2  molecules. 
From  the  results  of  these  PIMC  simulations  we  can  determine  the  average  energy,  as  a  function 
of  the  Al-Al  distance.  Since  the  Al  atoms  are  ~13  times  heavier  than  a  p\\2  molecule,  this 
average  energy  will  provide,  in  a  Born-Oppenheimer  sense,  the  effective  potential  energy  for 
motion  of  the  Al  atoms  along  the  recombination  coordinate.  If  a  potential  barrier  exists,  we  can 
then  use  standard  transition  state  theory  to  estimate  the  lifetime  against  recombination.  This 
approach  allows  the  investigation  of  atomic  recombination  in  solids  when  quantum  effects  are 
substantial. 
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For  substitution  sites  within  a  distance  of  ~  13  bohr,  we  found  that  the  two  A1  atoms 
greatly  distort  the  lattice  structure  to  allow  recombination  into  the  Al2  molecule.44  This  releases 
a  large  amount  of  energy.  This  result  is  different  from  the  model  studies  presented  in  a  previous 
paper,  where  the  pH2  molecules  were  fixed  in  place.  It  is  the  quantum  nature  of  the  pW2 
molecules  that  facilitates  recombination  of  the  A1  atoms.  However,  when  the  two  A1  atoms  are 
embedded  initially  at  distances  longer  than  ~14  bohr,  the  equilibrated  separation  between  the  two 
A1  atoms  is  little  changed  from  the  initial  substitutional  sites.  The  interaction  between  the  two 
A1  atoms  is  small,  and  the  A1  atoms  affect  the  system  separately. 

This  observation  from  our  simulations  provides  a  bound  to  the  maximum  concentration  to 
which  A1  atoms  could  be  doped  into  solid  p\\2,  without  recombination.  We  estimate  that 
recombination  will  occur  for  molar  ratios  greater  than  1  A1  :  15  H2,  or  0.0667  mole  percent. 

From  a  straightforward  application  of  simple  transition  state  theory,  we  estimated  the  lifetime  of 
the  embedded  solids  at  T  =  4  K  to  range  from  minutes  to  days.44 

G.  Investigation  of  Non-adiabaticity  in  the  F  +  HC1  — *  HF  +  Cl  and  0(3P )  +  H2  — *  OH  +  H 
Reactions 

Millard  Alexander  supervised  Michael  Deskevich,  a  graduate  student  at  the  University  of 
Colorado,  in  the  determination  of  the  coupled  electronic  potential  energy  surfaces  for  the  F  + 

HC1  — »  HF  +  Cl  reaction.  This  reaction  is  under  ongoing  experimental  investigation  in  the  group 
of  David  Nesbitt.  With  the  goal  of  future  scattering  calculations  to  compare  with  these 
experiments,  it  is  necessary  to  transform  the  three  electronically  adiabatic  potential  energy 
surfaces  (two  of  A'  symmetry  and  one  of  A"  symmetry)  into  a  quasi-diabatic  representation. 
This  is  entirely  analogous  to  our  earlier  treatment  of  the  F  +  H2  reaction  dynamics.45  Our 
expertise  in  the  determination  of  the  potential  energy  surfaces  for  the  F  +  H2  reaction,  as  well  as 
for  the  0(3/>)  H2  complex46  (the  latter  projected  completed  with  the  support  of  an  earlier 
AFOSR  grant),  proved  invaluable  in  mentoring  Mr.  Deskevich. 

We  also  extended  this  past  project  to  the  determination  of  the  full,  coupled  potential 
energy  surfaces  for  the  O (  V)  +  H2  — *  OH  +  H  reaction.  Here,  there  are  two  potential  energy 
surfaces  of  A"  symmetry  and  one  of  A'  symmetry.  In  contrast  to  the  F  +  H2  reaction,  two  of 
these  surfaces  lead,  over  a  high  barrier,  to  reaction,  while  one  is  repulsive.  Schatz  and  co¬ 
workers47  had  shown  that  the  crossing  with  the  singlet  potential  energy  surfaces  which  emanate 
from  the  higher  O(’D)  +  H2  occur  beyond  the  barrier  on  the  triplet  potential  energy  surfaces.  It 
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is  our  goal  to  carry  out  a  fully-quantum,  multiple-potential-energy  surface  study  of  the  O (3P)  + 
H2  —*  OH  +  H  reaction,  similar  to  our  earlier  studies  of  the  F  +  H2  and  Cl  +  H2  reactions.45,48 
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